Analysis of the spectral function of Ndi.gsCeo.isCuC^, obtained by angle resolved 

photoemission spectroscopy 
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Samples of Nd2- a; Ce a; Cu04, an electron-doped high temperature superconducting cuprate 
(HTSC), near optimal doping at x = 0.155 were measured via angle resolved photoemission 
(ARPES). We report a renormalization feature in the self energy ("kink") in the band disper- 
sion at ~ 50 — 60 meV present in nodal and antinodal cuts across the Fermi surface. Specifically, 
while the kink had been seen in the antinodal region, it is now observed also in the nodal region, 
reminiscent of what has been observed in hole-doped cuprates. 



.2,31 



The high temperature superconducting cuprates 
(HTSCs) have received much attention since their dis- 
covery, as the mechanism of their unusually high 
critical temperatures remains yet to be determined. 
Nd2- x Ce x Cu04 (NCCO) is a class of cuprates that re- 
sides on the relatively less-studied electron (n)-doped side 
of the phase diagram^, which is qualitatively different 
from the hole (p)-doped side. NCCO has recently been 
attracting increased interes^i 3 -^^^, in particular re- 
garding its links to their more commonly studied hole- 
doped counterparts^. 

Angle-resolved photoemission spectroscopy (ARPES) 
is used as a powerful direct probe of the electronic band 
structure (the one-particle spectral function)^. Sudden 
changes in the slope in the band dispersion observed via 
ARPES allow one to infer the underlying physics caus- 
ing these renormalizations^. Low-energy renormaliza- 
tions have been prevalently found in the nodal region of 
p-doped cuprates via ARPES 10 i n , but not in the nodal 
region of n-doped cuprates^, casting doubt on the uni- 
versality of the renormalization effect in the cuprates for 
the whole phase diagram. 

In this communication, we confirm a renormalization 
in NCCO near optimal doping (x = 0.155) around 
55 meV in the antinodal (X-M) region and report one 
in the nodal (r~M) region of about the same energy 
which was not observed earlier—. We do not observe 
any change of these renormalizations across the super- 
conducting (SC) phase transition. These results suggest 
that an oxygen phonon mode with comparable energy 
scale is a likely origin of these renormalizations. In light 
of the presence of phononic coupling effects, one needs to 
carefully exam the data to separate the respective contri- 
butions to the low energy spectra from AF band folding 
and lattice coupling effects, as the latter can also break 
the low energy dispersion into two "branches" . 

ARPES data were taken at beamline 5-4 of the Stan- 
ford Synchrotron Radiation Laboratory with Scienta 
SES200 and R4000 analyzers at a photon energy of 
16.75 eV. The a-(b-)axis was aligned at 45° to the light 
polarization for the nodal cut, and parallel to the light 
polarization for the other cuts. The energy resolution was 



« 10 meV and the angular resolution 0.3°. All sam- 
ples were cleaved at pressures better than 3 x 10~ n torr; 
the measurement temperature was 10 K unless other- 
wise specified. Single crystals of NCCO were grown at 
Stanford University^ 2 -. The doping level, determined by 
inductively coupled plasma spectroscopy, is 15.5%±0.7% 
Ce. Tc was determined by SQUID magnetometry to be 
25 K with a transition width of ss2 K. 

We show three different cuts through the Fermi surface 
(FS) in fig. [T] one in the antinodal, one in the nodal, and 
one near the hot spot region (region of low intensity at 
Ep near the crossing of the FS with the antiferromag- 
nctic (AF) Brillouin zone (BZ) boundary). The insets 
in fig. [1] a, c, and e show the respective locations of 
the three cuts in the BZ relative to the FS of NCCO. 
The raw data of these three cuts, including their EDC 
curves, are shown in fig. Q] a-f. MDC analysis, i.e. a 
fit with a Lorentzian lineshape, is used to extract Re £ 
and Im S: Im £ corresponds to the extracted width of 
the Lorentzian times the Fermi velocity, while Re £ is 
obtained by subtracting the extracted position from an 
assumed bare band. Apart from the assumptions and 
details of this analysi o 5 ' 13 , we phenomenologically apply 
MDC/EDC spectral function analysis as a tool for quan- 
tification, lacking a definitive theory. 

The raw spectrum in fig. Q] a of the antinodal cut re- 
veals a renormalization at about 50 meV binding energy. 
The corresponding EDC curves in fig. [T]b display a sharp 
quasiparticle peak near the Ep crossing. This peak ter- 
minates with a dip at about 50 meV binding energy The 
linewidth decreases upon approaching Ep. To determine 
the self energy, a parabolic barcband has been assumed 
and constructed by fixing the iJp-crossing point and the 
band bottom, which are determined by a linear fit to a 
MDC analysis < 40 meV, and a parabolic fit to an EDC 
analysis > 250 meV, respectively (cf. fig. [1] a). The 
band bottom thus obtained is at about 300 meV. The 
discrepancy between the EDC and MDC fitting around 
250 meV shows the limits of both analyses^. 

The raw spectrum of the nodal region cut depicted in 
fig. [T] c exhibits a similar renormalization feature at a 
similar energy of about 50 meV as seen in the antinodal 
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FIG. 1: (Color online) FS cuts and EDCs of the antinodal 
(a, b), nodal (c, d), and hot spot (e, f) regions. The insets 
in a), c), and e) depict a map of the upper right quadrant of 
the FS; the cut positions are indicated by pink (gray) lines, 
a), c), and e) show the raw spectra, while b), d), and f) show 
EDCs near the Fermi-crossing of the respective raw spectra 
on the left. Also indicated in the raw spectra are the MDC 
fits in solid black, EDC fits in dashed black, and assumed 
bare band in dotted red (black). The EDCs in b), d), and 
f) at Ef crossing are highlighted in thick red (black). A 
renormalization is seen in the nodal and antinodal regions 
around 50 meV (gray dashed arrows). 



region. The EDCs (fig. [T]d) show a peak at the Ep cross- 
ing. Again, the peak terminates at roughly 50 meV. A 
resemblance is found in the peak-dip-hump (PDH) struc- 
ture seen over a wide temperature range in the nodal 
EDCs of the single-layer compound Bi220lAi which is at- 
tributed to electron-phonon interaction, although a hump 
is absent here in the antinodal and faint in the nodal cuts. 
The assumed linear bare band was modeled by connect- 
ing the Ep crossing point obtained by a linear fit to the 
MDC dispersion < 40 meV, and a point at 300 meV 
obtained by a linear fit to the MDC dispersion between 
250 and 330 meV (cf. fig. Q]c). We note that assuming 
different bare bands changes the extracted renormaliza- 
tion strength, but not appreciably the energy of sharp 



features in the real part of the self energy. 

The spectrum of the hot spot region is displayed in fig. 
[T]e. Since the hot spot is located on the AF BZ bound- 
ary, the low-energy physics and thus the kink will be 
influenced by a crossover of two bands that could orig- 
inate from backfolding at the AF BZ boundary as we 
discuss later. Bare band assumption and construction 
are analogous to the antinodal region. 

The real and imaginary parts of the self energy are dis- 
played in fig. [2] a-f. A peak in Re £ is clearly visible at 
60 ± 5 meV in the antinodal region (fig. [2] a), with a cor- 
responding drop at about the same energy in Im £ (fig. 
Wp)i as expected from causality arguments. Likewise, in 
the nodal region a plateau is observed at about 50 ± 10 
meV in Re £ (fig. d]c). The drop in Im £ (fig. H]d) is 
much less clear, but the data seem qualitatively consis- 
tent with the Kramcrs-Kronig relations. The difference 
in our ability to extract the real and imaginary parts of 
the self-energy can easily be understood from the fact 
that, for a relatively broad feature, it is much easier to 
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FIG. 2: (Color online) Self energies extracted from the antin- 
odal (a, b), nodal (c, d), near-hot spot (e, f) cuts shown in fig. 
[H respectively, a), c), and e) show the respective real part; 
b), d), and f) the respective imaginary part. Changes in Re 
E and their corresponding drops in Im E are marked by green 
(gray) arrows. The error bars at selected points represent the 
3a confidence levels from the MDC-Lorentzian fit position in 
the Re E graphs, and the MDC-Lorentzian fit HWHM in the 
Im E graphs. The dashed arrows mark weak changes in Re E 
that are hard to distinguish from the statistical background. 
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FIG. 3: (Color online) Temperature dependence of the MDC- 
derived dispersions in the antinodal (a) and nodal (b) regions, 
with 3<r errors from the fit (thick horizontal bars). The black 
dotted lines serve as guide to the eye. The data have been 
average-subtracted in the gray region (dotted arrow) to over- 
lay them, c) Temperature dependent velocities in both nodal 
and antinodal region above and below the kink energy, ex- 
tracted from gray regions indicated by I-IV in a) and b). d) 
Temperature dependent fractions of the respective velocities 
below and above the kink energy for the nodal and antinodal 
region. Tc is indicated in c) and d) (vertical gray line). 



determine its position than width. 

We believe the reason why we are now able to resolve 
the kink in the nodal region and not before^ is mainly 
due to the improvement of the momentum resolution. 
The photon energy (16.75 eV) we used as compared to 
that used in previous measurements (53 eV) doubles our 
resolution in k-space. Indeed, the line width — which 
depends significantly on sample quality and resolution 
— is smaller (cf. fig. Ob, d) in our data presented in this 
communication. 

A change in slope can also be seen in Re £ near the 
hot spot region (fig. f5]e) around 50 ± 10 meV. However, 
the relation of its cause to the features in the nodal and 
antinodal regions remains speculative, even though we 
find its appearance at the same energy suggestive. 

Lastly, a weak change of slope at a higher energy of 
around 160 meV can be observed in Re £ in both the 
nodal and antinodal region (fig. [2]a,c). Within the statis- 
tics of the data, no change is visible in either of the corre- 
sponding Im £ (fig. f2]b,d) at this energy. It is currently 
being investigated and will not be discussed here. 

The temperature dependence of the observed ubiqui- 
tous feature in the nodal and the antinodal region is 
shown in fig. 3. Within the statistics, there is no change 



visible in cither region below (<6 K, 10 K) and above 
(30 K, 50 K) the SC transition. 

For further clarification, the velocities obtained by 
linear fits against the MDC Lorentz position above (- 
40 meV< E - E F <-20 meV) and below (-135 meV< 
E — Ep <-100 meV) the renormalization are graphed 
versus temperature for both nodal and antinodal regions 
in fig. [3] c. The data result from seven different mea- 
surements on four samples in the nodal, and four differ- 
ent measurements on three samples in the antinodal re- 
gion. First, neither of these velocities change appreciably 
across the SC transition. Second, the Fermi velocities of 
the nodal and antinodal regions are the same, confirming 
previous results^. Third, the ratio of renormalized and 
unrenormalized velocities — which are indicative of the 
renormalization strength of the electronic band — do not 
seem to change within the statistics in both nodal and 
antinodal regions, cf. fig. [3]d. 

The n-doped HTSCs possess hot spots that lie on the 
AF BZ boundary^, and recently, pockets were observed 
suggesting a band back-folding at the AF BZ boundary^; 
a possible mechanism could be (7r, 7r)-scattering due to 
some short range or remnant antiferromagnetic (SDW) 
ordering. To first order, this folding results in a pair of 
new band sheets (cf. fig. 0] a)^. Whenever there is a 
crossover from one folded band sheet to another, as for 
example at (7r, 0) in the antinodal region (solid curve to 
solid curve in fig. [5]), one expects a kink- like feature in the 
dispersion. We do not however think that this is the sole 
cause for the low energy feature, since this crossover takes 
place at the band bottom (w 300 meV) in the antinodal 
and above the Fermi energy in the nodal region, leading 
to a highly anisotropic behavior both in shape and en- 
ergy of the feature, contradicting our observed isotropy 
in both energy and renormalization strength. 

If we were to invoke an electron energy dependence 
of the scattering, V^^i^h), whereas the interaction is 
nonzero only for \ujk\ < u max ps 50 meV, a fade-over 
is produced between the unrenormalized, unfolded band 
(dotted curve in fig. [4j and one of the renormalized, 
folded bands (solid curves; cf. fig. 2] b and c) around 
Umax- However, we do not know of any interaction with 
above-mentioned energy dependence. Provided it exists, 
the produced crossover results in a downwards kink in 
the nodal and an upwards kink in the antinodal region, 
contradicting our observations. 

We also argue against a magnetic resonance or spin 
flip waves as the origin of the low energy renormaliza- 
tion we discuss here, since the characteristic energy scales 
for electron-doped cupratcs are too different from the 
55 meV observed by us. First, the relevant region of 
the spin excitation dispersion relation is about 0.3 eV— . 
Second, a resonance has been reported at about 10 meV 
in PLCCOil and in NCCOH. This is still disputed; 
at least for NCCO, recent data show even lower energy 
scales (within 4-8 meV) for the SC gap and a possible 
resonance^. Whether a more complicated interaction 
mechanism could explain our data is to be investigated. 
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FIG. 4: (Color online) Schematic tight binding (TB) models. 
Solid blue (black) curves correspond to TB without, dotted 
red (black) ones to TB with simple perturbative V^^ scat- 
tering, a) cuts along high symmetry directions in the BZ. b) 
antinodal and c) nodal details of a). The enlarged regions 
shown in b) and c) are indicated by rectangles in a). 



In addition, the fact that the kink at both nodal and 
antinodal region does not change across the SC transi- 
tion lends further support to the electron-phonon cou- 
pling scenario. The absence of a change of the kink en- 
ergy across the superconducting phase transition is due 
to the small superconducting gap in NCCO. Also, the 
EDCs at the Ep crossing exhibit the same structure as 
in Bi2201 — a small peak near Ep and a shallow dip at 
the same energy of the rcnormalization seen in Re £ (cf. 
fig.fflb,d)U. * 

Within the context of electron-phonon interaction, a 
rcnormalization at around 70 meV is seen in the nodal 
region of p-dopcd HTSCs in ARPES and has been at- 
tributed to electron-phonon interaction o 10 i 26 . A soften- 
ing of an oxygen mode was seen as well for the n-doped 
compound NCCO via inelastic neutron and X-ray scat- 
tering (55 meV) 23 ' 27 , and now also in the nodal region 
via ARPES as observed by us in the same relevant en- 
ergy region. Independently, Park et alj2£ and Liu et al^i 
discovered a kink in the nodal region in several n-doped 
HTSCs and came to very similar conclusions regarding its 
origin^. If interpreted as due to electron-phonon cou- 
pling, our results suggest that electronic coupling to the 
oxygen phonon mode may be even more universal, now 
stretching across the phase diagram. We speculate that 
these phonon modes may play a vital role in determining 
the low energy physics of the HTSC o 20 ' 22 . 



A phononic origin on the other hand could account 
for the isotropy of the renormalization strength and en- 
ergy of the kink observed by us in a simple picture (al- 
though more intricate phononic interactions may result 
in anisotropic renormalization o 20 ' 21 and even d-wave SC 
gap a 20 ' 22 ). Furthermore, a softening with doping of oxy- 
gen phonon modes in the relevant energy regime has been 
observed via inelastic x-ray scattering^ and inelastic 
neutron scattering 2 ^. Raman and IR spectra also show 
the existence of oxygen modes of the relevant energies 2 ^. 
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